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A classical model is presented for the features of parametric down-conversion and homodyne 
detection relevant to recent proposed "loophole-free" Bell tests. The Bell tests themselves are 
uncontroversial: there are no obvious loopholes that might cause bias and hence, if the world does, 
after all, obey local realism, no violation of a Bell inequality will be observed. Interest centres 
around the question of whether or not the proposed criterion for "non-classical" light is valid. If 
cit is not, then the experiments will fail in their initial concept, since both quantum theorists and 
local realists will agree that we are seeing a purely classical effect. The Bell test, though, is not 
the only criterion by which the quantum-mechanical and local realist models can be judged. If the 
experiments are extended by including a range of parameter values and by analysing, in addition 
to the proposed digitised voltage differences, the raw voltages, the models can be compared in their 
overall performance and plausibility. 
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I. INTRODUCTION 

No test of Bell's inequalities 0,01 to date has been free 
of "loopholes". This means that, despite the high levels 
of statistical significance frequently achieved, violations 
of the inequalities could be the effects of experimental 
bias of one kind or another, not evidence for the pres- 
ence of quantum entanglement. Recent proposed exper- 
iments by Garci'a-Patron Sanchez Q et al. and Nha and 
Carmichael ^ show promise of being genuinely free from 
such problems. If the world in fact obeys local realism, 
they should not, therefore, infringe any Bell inequality. 




FIG. 1: Proposed experimental set-up. In the current text, 
phase shifts 9 and cj) are renamed 9a and 9b. {Reproduced 
with permission from Sanchez et al., Phys. Rev. Lett. 93, 
130409 (2004)) 



'Electronic address: 'ch.th ompsonl@virgin.iietl 



The current article presents a classical (local realist) 
model that should be able, once all relevant details are 
known, to explain the results. It uses standard classical 
theory for homodyne detection, but for the behaviour 
of the "Optical Parametric Amplification" ("Parametric 
Down-Conversion") source introduces new ideas. 

As far as the "loophole-free" status of the proposed 
experiments is concerned, there would appear to be no 
problem. A difficulty that seems likely to arise, though, 
is that theorists may not agree that the test beams used 
were in fact "non-classical", so the failure to infringe a 
Bell inequality will not in itself be interpreted as show- 
ing a failure of quantum mechanics 5]. The criterion 
to be used to establish the non-classical nature of the 
light is the observation of negative values of the Wigner 
density, and there is reason to think that, even if the 
standard method of estimation seems to show that these 
are achieved, this may not in fact be so. 

In the quantum mechanical theory discussed in the 
proposals and in other recent papers ]^ Q , the presence 
of negative Wigner densities is seen almost as a corollary 
of the observation of double peaks in the distribution of 
the phase-averaged voltage differences recorded by the 
homodyne detectors. Such double peaks, however, can 
readily be shown to occur naturally, under purely classi- 
cal assumptions, so long as noise levels are low. They are 
a consequence of the fact that a sine function spends pro- 
portionately more time near the extremes than around 
the central, zero, value. It would seem that the theory 
that leads from their observation to the deduction of non- 
classicality may be in error. This, though, is a problem 
that is of no direct relevance to the classical model under 
consideration. Far from being, as suggested by Sanchez 
and others, the "hidden variable" needed, Wigner density 
plays no part whatsoever. 

Regardless of the outcome of the Bell tests, and 
whether or not the light is declared to be non-classical, 
there are features of the experiments that can usefully 
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be exploited to compare the strengths of quantum me- 
chanics versus (updated) classical theory as viable mod- 
els. The two theories approach the situation from very 
different angles. Classical theory traces the causal re- 
lationships between phenomena, starting with the sim- 
plest assumption and building in random factors later 
where necessary. Quantum mechanics starts with mod- 
els of complete ensembles, all random factors included. 
This, it is argued, is inappropriate, since two features of 
the proposed experiments demand that we consider the 
behaviour of individual events, not whole ensembles: the 
process of homodyne detection itself, and the Bell test. 



II. THE PROPOSED EXPERIMENTS 

The experimental set-up proposed by Sanchez et al. 
is shown in Fig. 1, that of Nha and Carmichael being 
similar. In the words of the Sanchez et al. proposal: 

The source (controlled by Sophie) is based 
on a master laser beam, which is converted 
into second harmonic in a nonlinear crystal 
(SHG). After spectral filtering (F), the sec- 
ond harmonic beam pumps an optical para- 
metric amplifier (OPA) which generates two- 
mode squeezed vacuum in modes A and B. 
Single photons are conditionally subtracted 
from modes A and B with the use of the beam 
splitters BS^ and BS^ and single-photon de- 
tectors PD^ and PDg. Alice (Bob) mea- 
sures a quadrature of mode A (B) using a 
balanced homodyne detector that consists of 
a balanced beam splitter BS3 (BS4) and a 
pair of highly-efficient photodiodes. The lo- 
cal oscillators LO^i and LOb are extracted 
from the laser beam by means of two addi- 
tional beam splitters BSi and BS2. 

The classical description, working from the same fig- 
ure, is just a little different. Quantum theoretical terms 
such as "squeezed vacuum" and are not used since they 
are not appropriate to the model and would cause con- 
fusion. The description might run as follows: 

The master laser beam (which is, incidentally, pulsed) 
is frequency-doubled in the crystal SHG. After filtering 
to remove the original frequency, the beam is used to 
pump the crystal OPA, which outputs pairs of classical 
wave pulses at half the input frequency, i.e. at the original 
laser frequency. The selection of pairs for analysis is done 
by splitting each output at an unbalanced beamsplitter 
(BSa or BSb), the smaller parts going to sensitive de- 
tectors PDa or PDb- Only if there are detections at 
both PDa and PDs is the corresponding homodyne de- 
tection included in the analysis. The larger parts proceed 
to balanced homodyne detectors, i.e. ones in which the 
intensities of local oscillator and test inputs are approxi- 
mately equal. The source of the local oscillators LOa and 



LOb is the same laser that stimulated, after frequency 
doubling, the production of the test beams. 



III. HOMODYNE DETECTION 

In (balanced) homodyne detection, the test beam is 
mixed at a beamsplitter with a local oscillator beam of 
the same frequency and the two outputs sent to pho- 
todetectors that produce voltage readings for every input 
pulse. In the proposed "loophole-free" Bell test the dif- 
ference between the two voltages will be converted into 
a digital signal by counting all positive values as -1-1, all 
negative as -1. 
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FIG. 2: Inputs and outputs at the beamsplitter in a homo- 
dyne detector. El is the local oscillator beam, E the test 
beam, Et and Er the transmitted and reflected beams respec- 
tively. 

Assuming the inputs are all classical waves of the same 
frequency and there are no losses, it can be shown (see 
below) that the difference between the intensities of the 
two output beams is proportional to the product of the 
two input intensities multiplied by sin^, where 9 is the 
phase difference between the test beam and local oscil- 
lator. If voltages are proportional to intensities then it 
follows that the voltage difference will be proportional 
to sin 9. When digitised, this transforms to a step fimc- 
tion, taking the value —1 for — tt < 9 < and +1 for 
< 9 < TT. (The function is not well defined for integral 
multiples of tt.) 



A. Classical derivation of the predicted voltage 
difference 

Assume the test and local oscillator signals have the 
same frequency, w, the time-dependent part of the test 
signal being modelled by e"^ , where (ignoring a constant 
phase offset 8]) (f> = ojt is the phase angle, and the local 
oscillator phase and test beam phases differ by 9. [Note 
that although complex notation is used here, only the 
real part has meaning: this is an ordinary wave equation, 
not a quantum-mechanical "wave function" . To allay any 
doubts on this score, the derivation is partially repeated 
with no complex notation in the Appendix.] 
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Let the electric fields of the test signal, local oscillator 
and reflected and transmitted signals from the beamsplit- 
ter have amplitudes E, El, E^ and Et respectively, as 
shown in Fig. [21 Then, after allowance for phase delays 
of 7r/2 at each reflection and assuming no losses, we have 



v2 



and 



Elc 



(1) 



(2) 



The intensity of the reflected beam is therefore 
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Similarly, it can be shown that the intensity of the trans- 
mitted beam is 



EtE; = ^{E^ +El- 2EEl sin( 



(4) 



If the voltages registered by the photodetectors are pro- 
portional to the intensities, it follows that the difference 
in voltage is proportional to 2EEl sin 6*. When digitised, 
this translates to the step function mentioned above. The 
probabilities for the two possible outcomes are, as shown 
in Fig. 13 



and 
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1 for -TT < 61 < 
for < 61 < vr 



for -TT < 6* < 

1 for < 61 < vr 



(5a) 



(5b) 



Note that the probabilities are undefined for integral 
multiples of vr. In practice it would be reasonable to 
assume that, due to the presence of noise, all the values 
were 0.5, but for the present purposes the integral values 
will simply be ignored. 



IV. FRESH IDEAS ON PARAMETRIC 
DOWN-CONVERSION 

If the frequencies and phases of both test beams and 
both local oscillators were all identical apart from the 
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FIG. 3: Probabilities of and '-' outcomes versus phase 
diflerence, using digitised difference voltages from a perfect, 
noise-free, balanced homodyne detector. 



applied phase shifts, the experiment would be expected to 
produce step function relationships between counts and 
applied shifts both for the individual (singles) counts and 
for the coincidences. 

It may safely be assumed that this is not what is ob- 
served. It would have shown up in the preliminary trials 
on the singles counts (see ref. Q), which would have fol- 
lowed something suggestive of the basic predicted step 
function as the local oscillator phase shift was varied. 
What is observed in practice is more likely to be simi- 
lar to the results obtained by Schiller et al. Their 
Fig. 2a, reproduced here as Fig. 0] shows a distribution 
of photocurrcnts that is clustered around zero, for 9 tak- 
ing integer multiples of tt, but is scattered fairly equally 
among positive and negative values in between. 
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FIG. 4: A typical scatter of "noise current" (related to volt- 
age difference) for varying local oscillator phase. {Reproduced 
from S. Schiller et al., Phys. Rev. Lett. 87 (U), 2933 (1996). 
Permission requested, June 12, 2005.) 

When digitised, the distribution would reduce to two 
straight horizontal lines, showing that for each choice of 9 
there is an equal chance of a or a '— ' result. As in any 
other Bell test setup, though, the absence of variations in 
the singles counts does not necessarily mean there is no 
variation in the coincidence rates. However, as explained 
in the next section, the predicted coincidence curves are 
not the zig-zag ones of standard classical theory. These 
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would be expected if we had full "rotational invariance" 
[IQi] . If the ideas presented here are correct, we have in- 
stead, in the language of an article by the present author 
[ll], only binary rotational invariance. Schiller's scatter 
of photocurrent differences is seen as evidence that the 
relative phase can (under perfect conditions) take just 
one of two values, or tt. The scatter is formed from 
a superposition of two sets of points, corresponding to 
two sine curves that are out of phase, together with a 
considerable amount of noise. 

It is suggested that the two "phase sets" arise from the 
way in which the pulses are produced, which involves, af- 
ter the frequency doubling, the degenerate case of para- 
metric down-conversion, the latter producing pulses that 
are (in contrast to the quantum-mechanical assumption 
of conjugate frequencies) of exactly equal frequency. Con- 
sider an initial pump laser of frequency uj. In the pro- 
posed experiment, this will be doubled in the crystal SHG 
to 2llj then halved in OPA back to uj. At the frequency 
doubling stage, one laser input wave peak gives rise to 
two output ones. Assuming that there are causal mecha- 
nisms involved, it seems inevitable that every other wave 
peak of the output will be exactly in phase with the in- 
put. When we now use this input to produce a down- 
conversion, the outputs will be in phase either with the 
even or with the odd peaks of the input, which will make 
them either in phase or exactly out of phase with the 
original laser. 

We thus have two classes of output, differing in phase 
by TT. If we define the random variable a to be for one 
class, TT for the other, this will clearly be an important 
"hidden variable" of the system. 

The theory of (degenerate) parametric down- 
conversion assumed here differs in several respects from 
the accepted quantum-mechanical one, and also from 
Stochastic Electrodynamics. No attempt is made to give 
a full explanation of the physics involved. The theory 
is rather of the nature of an empirical result, effectively 
forced on us by a number of different experiments 
in quantum optics (to be discussed in later papers). 
Accepted theory says (if, indeed, it allows at all for 
the existence of "photons" with definite frequencies 
p^) that though the sum of the frequencies of the two 
outputs is equal to that of the pump laser, even in the 
degenerate case the two will differ slightly 13]. The 
exact frequency of one will be |(wo + Sur), that of the 
other ^{ujq ~ where Wq is the pump frequency. If 
this is the case in the proposed experiment, though, 
it will severely reduce the visibility of any coincidence 
curve observed when the experimental beam is mixed 
back with the source laser in the homodyne detector. 

The preliminary experiments (see ref. [^) using just 
one output beam may already be sufficient to show that 
the interference is stronger than would be the case if there 
were any difference between local oscillator and test beam 
frequencies. It is known that the source laser has quite a 
broad band width, i.e. that luq is not constant. Though 
it is likely that it is only part of the pump spectrum that 



induces a down-conversion, so that the band width of 
the test beam may be considerably narrower than that of 
the pump, it too is non zero. It follows that agreement 
of frequency between this and the test beam must be 
because we are always dealing, in the degenerate case, 
with exact frequency doubling and halving. 

V. A CLASSICAL MODEL OF THE PROPOSED 
BELL TEST 

In the proposed Bell test of Sanchez et al., positive 
voltage differences will be treated as -1-1, negative as - 
1. Applying this version of homodyne detection to both 
branches of the experiment, the CHSH test (—2 < 5 < 2) 
will then be applied to the coincidence counts. Under 
quantum theory it is expected that, so long as "non- 
classical" beams are employed, the test will be violated. 
However, since there are no obvious loopholes in the ac- 
tual Bell test (see later), local realism should win: the 
test should not be violated. In the classical view, this 
prediction is unrelated to any supposed non-classical na- 
ture of the light. 

A. The basic local realist model 

If we take the simplest possible case, in which to all 
intents and purposes all the frequencies involved are the 
same, the hidden variable in the local realist model is 
clearly going to be the phase difference (a = or tt) 
between the test signal and the local oscillator. If high 
visibility coincidence curves are seen, it must be because 
the values of a are identical for the A and B beams. As- 
suming no noise, the basic model is easily written down. 

From equation (|5a|l . the probability of a —1 outcome 
on side A is 

fa \ / 1 for -TT < - a < 
p.(0^,a) = |Q forO<0^-a<7r, 

where 6a is the phase shift applied to the local oscillator 
A, a is the hidden variable and all angles are reduced 
modulo 2tt. Similarly, the probability of a -1-1 outcome is 

, . / for -TT < 6*^ - a < , . 

p+(0^,a) = |^ forO<0^-a<^, 

Assuming equal probability ^ for each of the two possi- 
ble values of a JJl , the standard "local realist" assump- 
tion that independent probabilities can be multiplied to 
give coincidence ones leads to a predicted coincidence 
rate of 

+ ^p+{9A,Tr)p+{9B,7T), (8) 
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with similar expressions for 



and P__. 



The resuh for 6^ — 7r/2, for example, is 



P++(7r/2,0B) = 
For Oa = -tt/2 it is 



for -TT kOb <0 
1/2 for 0<9b <n. 



/ ' -"^ I ior < 9b < n. 



(9) 



(10) 



Note that, as illustrated in Fig.|Sl the coincidence prob- 
abilities cannot, in this basic model, be expressed as func- 
tions of the difference in detector settings, 9b — 9 a- This 
failure, marking a significant deviation from the quantum 
mechanical prediction, is an inevitable consequence of the 
fact that we have (as mentioned earlier) only binary, not 
full, rotational invariance. 

(a)P++(,:/2,eB) 
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(b)P++(_;:/2,eB) 
0.5 
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(c)P++(7i/4,eB) 
0.5 



0.0 



FIG. 5: Predicted coincidence curves for the ideal experiment, 
(a) and (b) illustrate the settings most likely to be chosen in 
practice, giving the strongest correlations. 9a is fixed at 7r/2 
or — 7r/2 while 9b varies. In theory, any value of 9a between 
and TT would give the same curve as (a), any between — tt 
and the same as (b). An example is shown in (c), where 
9a is 7r/4 but the curve is identical to (a). We do not have 
rotational invariance: the curve is not a function of 9b — 9a- 



B. Fine-tuning the model 

Many practical considerations mean that the final lo- 
cal realist prediction will probably not look much like the 
above step function. It may not even be quite periodic. 
The main logical difference is that, despite all that has 
been said so far, the actual variable that is set for the lo- 
cal oscillators is not directly the phase shift but the path 
length, and, since the frequency is likely to vary slightly 
from one signal to the next (though always keeping the 
same as that of the pump laser) , the actual phase differ- 
ence between test and local oscillator beams will depend 



on the path length difference and on the frequency. In 
a complete model, therefore, the important parameters 
will be path length and frequency, with phase derived 
from these. 

If frequency variations are sufhciently large, the sit- 
uation may approach one of rotational invariance (RI), 
but it seems on the face of it unlikely that this can be 
achieved without loss of correlation. If we do have RI, 
the model becomes the standard realist one in which the 
predicted quantum correlation varies linearly with differ- 
ence in phase settings, but it is more likely that what 
will be found is curves that are not independent of the 
choice of individual phase setting. They will be basically 
the predicted step functions but converted to curves as 
the result of the addition of noise. 

P++(7t/2,eB) 

0.5 n 



0.0 



FIG. 6: Likely appearance of coincidence curves in a real 
experiment with moderate noise. 

It is essential to know the actual experimental condi- 
tions. Several relevant factors can be discovered by care- 
ful analysis of the variations in the raw voltages in each 
homodyne detection system. If noise is low, the pres- 
ence of the two phase sets, and whether or not they are 
equally represented, should become apparent. All this 
complexity, though, has no bearing on the basic fact of 
the existence of a hidden variable model and the con- 
sequent prediction that the CHSH Bell test will not be 
violated. 



C. The role of the "event-ready" detectors 

In the quantum-mechanical theory, the expectation of 
violation of the Bell test all hinges on the production 
of "non-classical" light. The light directly output from 
the crystal OPA is assumed to be Gaussian, i.e. it takes 
the form of pulses of light that have a Gaussian intensity 
profile and also, as a result of Fourier theory, a Gaus- 
sian spectrum. When this is passed through an unbal- 
anced beamsplitter (BSa or BSb) and a "single photon" 
detected by an "event-ready" detector, the theory says 
that the subtraction of one photon leaves the remaining 
beam "non-Gaussian" . Although there is mention here 
of single photons, the theory is concerned with the en- 
semble properties of the complete beams, not with the 
individual properties of its constituent pulses. 

In the local realist (classical) model, the shapes of the 
spectra are not relevant except insofar as a narrow band 
width is desirable for the demonstration of dramatic cor- 
relations. The event-ready detectors play, instead, the 
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important role of selecting for analysis only the strongest 
down-converted output pairs, it being assumed that the 
properties of the transmitted and reflected light at the 
unbalanced beamsplitters are identical apart from their 
amplitudes. It is likely that those detected signals that 
are coincident with each other will be genuine "degener- 
ate" ones, i.e. of exactly equal frequency, quasi-resonant 
with the pump laser. The unbalanced beamsplitters and 
the detectors PD^i and PD^ need to be set so that the 
intensity of the detected part is sufficient to be above the 
minimum for detection but low enough to ensure that all 
but the strongest pulses are ignored. 

In neither theory are the event-ready detectors really 
needed in their "Bell test" role of ensuring a fair test (see 
below), since the homodyne detectors are almost 100% 
efficient. 



VI. VALIDITY OF THE PROPOSED BELL TEST 

Coincidences between the digitised voltage differences 
will be used in the CHSH Bell test 0, 0| , but avoiding 
the "post-selection" that has, since Aspect's second ex- 
periment p^ . become customary. The Sanchez et al. 
proposal is to use "event-ready" detectors, as recom- 
mended by Bell himself for use in real experiments [l8j |. 
None of the usual loopholes are expected to be ap- 
plicable: 

1. With the use of the event-ready detectors, non- 
detections are of little concern. The detectors 
(PD^ and PD^ in Fig. act to define the sam- 
ple to be analysed, and the fact that they do so 
quite independently of whether or not any mem- 
ber of the sample is then also detected in coin- 
cidence at the homodyne detectors ensures that 
no bias is introduced here. The estimate of 
"quantum correlation" to be used in calcu- 
lating the CHSH test statistic is E ^ {N++ab + 
N__AB - N+_AB - N^+ab)/Nab, where the iV's 
are coincidence counts and the subscripts are self- 
explanatory. This contrasts with the usual method, 
in which the denominator used is not Nab but the 
sum of observed coincidences, N^+ab + N^^ab + 

AB+N \-AB- The use of the latter can readily 

be shown to introduce bias unless it can be assumed 
that the sample of detected pairs is a fair one. That 
such an assumption fails in some plausible local re- 
alist models has been known since 1970 or earlier 

MM- 

2. There is no possibility of synchronisation problems, 
since a pulsed source is used. 

3. No "accidentals" will be subtracted. 

4. The "locality" loophole can be closed by using long 
paths and a random system for choosing the "de- 
tector settings" (local oscillator phase shifts) during 
the propagation of the signals. 



The system is almost certainly not going to be "ro- 
tationally invariant" (not all phase differences will be 
equally likely), but this will not invalidate the Bell test. 
It may, however, be important in another way. It is likely 
that high visibilities will be observed in the coincidence 
curves (i.e. high values of (max - min) / (max + min) in 
plots of coincidence rate against difference in phase shift), 
leading to the impression that the Bell test ought to be 
violated. These visibilities, though, will depend on the 
absolute values of the individual settings. High ones will 
be balanced by low, with the net effect that violation 
does not in fact happen. 

VII. VALIDITY AND SIGNIFICANCE OF 
NEGATIVE ESTIMATES FOR WIGNER 
DENSITIES 

Part of the evidence that is put forward as indicating 
that negative Wigner densities are likely to be obtained 
consists in the observation that, when 9 is varied ran- 
domly, the distribution of observed voltage differences 
shows a double peak (see Fig. [7|). There is a tendency 
to observe roughly equal numbers of -I- and - results but 
relatively few near zero. The fact that the relationship 
depends on the sine of 9 is, however, sufficient to explain 
why this should be so. 
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FIG. 7: Observed distribution of voltage differences, X, using 
homodyne detection in an experiment similar to the proposed 
Bell test and averaging over a large number of applied phase 
shifts. {Based on Fig. J^a of A. I. Lvovsky et al, Phys. Rev. 
Lett, 87, 050402 (2001).) 

To illustrate, let us consider the following. The sine 
of an angle is between and 0.5 whenever the angle is 
between and 7r/6. It is between 0.5 and 1.0 when the 
angle is between 7r/6 and 7r/2. Since the second range of 
angles is twice the first yet the range of the values of the 
sine is the same, it follows that if all angles in the range 
to TT are selected equally often there will be twice as 
many sine values seen above 0.5 as below. 

When allowance is made for the addition of noise, the 
production of a distribution such as that of Fig. [7| for 
the average when angles are sampled uniformly comes 
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as no surprise. Clearly, as the experimenters themselves 
recognise, the dip is not in itself sufficient to prove the 
non-classical state of the light. For this, direct measure- 
ment of the Wigner density is required, but there is a 
problem here. No actual direct measurement is possible, 
so it has to be estimated, and the method proposed is 
the Radon transformation Q^. It is claimed [23| that 
in other experiments Wigner densities calculated by this 
procedure have shown negative regions, but perhaps the 
method should be checked for validity? 

In any event, as already explained, the natural hidden 
variable relevant to the proposed experiment is the phase 
of the individual pulse, not any statistical property of the 
whole ensemble. 



VIII. SUGGESTIONS FOR EXTENDING THE 
EXPERIMENT 

The basic set-up would seem to present an ideal op- 
portunity for investigation of some key aspects of the 
quantum and classical models, as well as the operation 
of the Bell test "detection loophole" . 

The operation of the detection loophole could 
be illustrated if, instead of using the digitised differ- 
ence voltages of the homodyne detectors, the two sep- 
arate voltages are passed through discriminators. The 
latter operate by applying a threshold voltage that can 
be set by the experimenter and counting those pulses 
that exceed it. These can be used in a conventional 
CHSH Bell test, i.e. using total observed coincidence 
count as denominator in the estimated quantum correla- 
tions E. The model that has been known since Pearle's 
time (1970) predicts that, as the threshold voltage used 
in the discriminators is increased and hence the number 
of registered events decreased (interpreted in quantum 
theory as the detection efficiency being decreased), the 
CHSH test statistic S, if calculated using estimates E = 
{N++ + N_.-N+_-N_+)/{N++ + N^.+N+_ + N_+), 
will increase. If noise levels are low, it may well exceed 
the Bell limit of 2. 

The existence of the two phase sets could also 
be investigated if either the raw voltages or the undigi- 
tised difference voltages are analysed. So long as the 
noise level is low, the existence of the two superposed 
curves, one for a = and the other for a — tt, should 
be apparent. It would be interesting to investigate how 
the pattern changed as optical path lengths were varied. 
Schiller's pattern might be hard to reproduce using long 
path lengths, where exact equality is needed unless the 
light is monochromatic. 

If the primary goal of the experimenter is clearly set 
out to be the comparison of the performance of the two 
rival models, rather than merely the conduct of a Bell 
test, further ideas for modifying the set-up will doubtless 
emerge when the first experiments have been done. 



IX. CONCLUSION 

The proposed experiments would, if the "non- 
classicality" of the light could be demonstrated satisfac- 
torily, provide a definite answer one way or the other re- 
garding the reality of quantum entanglement. They could 
usefully be extended to include empirical investigations 
into the operation of the Bell test detection loophole. 
Perhaps more importantly, though, they present valuable 
opportunities to compare the performance of the two the- 
ories in both their total predictive power and their com- 
prehensibility. Are parameters such as "Wigner density" 
and "degree of squeezing" really the relevant ones, or 
would we gain more insight into the situation by talking 
only of frequencies, phases and intensities? Parameters 
such as the detection efficiency and the transmittance 
of the beamsplitters will undoubtedly affect the results, 
but do they do this in the way the quantum theoretical 
model suggests? It will take considerably more than just 
the minimum runs needed for the Bell test if we are to 
find the answers. 

The detailed predictions of the local realist model can- 
not be given until the full facts of the experimental set-up 
and the performance of the various parts are known, but 
it gives, in any event, a simple explanation of the double- 
peaked nature of the distribution of voltage differences. 
The peaks arise naturally from the way in which homo- 
dyne detection works, and the quantum theoretical idea 
that they are one of the indications of a non-classical 
beam or of negative Wigner density would not appear 
to be justifiable. The idea that a classical beam can be- 
come non-classical by the act of "subtracting a photon" 
is, equally, of doubtful validity. In the classical model, 
the only effect of the subtraction and detection of part 
of each beam is to aid the selection for coincidence anal- 
ysis of those pulses that are likely to be most strongly 
correlated. 
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APPENDIX A: ALTERNATIVE CLASSICAL 
DERIVATION OF THE HOMODYNE 
DETECTION FORMULA 

A derivation is given here that does not involve com- 
plex numbers and hence confirms that the equations 
in the text are ordinary wave equations, not quantum- 
mechanical "wave functions" . 

The relationship between intensity difference and the 
local oscillator phase can be checked as follows: 

Assume the two input beams are 
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Experimental beam: E cos <p , where, as before, (j> = wt 

Local oscillator: El cos{^ + 6) 

Then the output beams, assuming a 50-50 beamsplitter 
and no losses, can be written 

Reflected beam: 

Er = ■^{Ecos{(l) + 'K/2) + ELCOs{(j) + e)) 
V2 

= -^{-Esm.(t) + ELCOs{(t) + e)) (Al) 
v2 



the homodyne detector is therefore expected to be zero. 
A zero difference is also found for = tt, but other values 
of 6 produce more interesting results. 

For example, for = 7r/2 we find 



Er = — z={E/ sin ip){— sin ip sin (/) — cos tp sin (p) 
v2 

= — ^(S/sin'^) sin(/)(— sin^ — cosV'), (A4) 
v2 



Transmitted beam: 



Et = -j={Ecos(j) + ELCOs{(t) + e + TT/2)) 



%/2 



{Ecoscj}- ELsin{(j) + e)). 



(A2) 



Let us define a (constant) angle ip such that tau'^ = 
E/El, making Ecostjj = El sintjj. 

Consider the case when ^ = 0. We have 



Er = — ^(S/sin^)(— sin^sin^ + cosV'cos^) 
v2 



V2 



[E/ sinij)) cos{tp + (p), 



(A3) 



so that amplitude is proportional to -^E/ sin^ and in- 
tensity to ^E^/ sin^ ip. 

Similarly, the intensity of the transmitted beam is also 
proportional to sin^ tp. The voltage difference from 



Et = —1= {E/ sin ip) (sin ip cos cp — cos tp cos cp) 
v2 

= — ^(ii^/sin^) cos^(sin^ — cos V')- (A5) 
v2 



The difference in intensities is therefore proportional to 
Difference = ^ {E"^/ sin^ ip) [(- sin tp - cos ip)'^ 



— (sin V' — cos V')^] 

= ^{E'^/sin'^ip)4:Sinipcosip) 

= —E"^ cos Ip / simp 
= 2E'^El/E 

= 2EEl. (A6) 

This is consistent with the result obtained by the 
method in the main text, so it seems safe to accept that 
as being correct. 
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